Preprint typeset in JHEP style - HYPER VERSION 



February 2, 2008 
MSU-HEP-xxxxxx 



Uncertainties of the Inclusive Higgs Production Cross 
Section at the Tevatron and the LHC 



Alexander Belyaev, Jon Pumplin, Wu-Ki Tung, C.-P. Yuan 

Department of Physics and Astronomy, Michigan State University, East Lansing, MI 
48824, USA 



Abstract: We study uncertainties of the predicted inclusive Higgs production cross section 
due to the uncertainties of parton distribution functions (PDF). Particular attention is 
given to bbH Yukawa coupling enhanced production mechanisms in beyond SM scenarios, 
such as MSSM. The PDF uncertainties are determined by the robust Lagrange Multiplier 
method within the CTEQ global analysis framework. We show that PDF uncertainties 
dominate over theoretical uncertainties of the perturbative calculation (usually estimated 
by the scale dependence of the calculated cross sections), except for low Higgs masses at 
LHC. Thus for the proper interpretation of any Higgs signal, and for better understanding 
of the underlying electroweak symmetry breaking mechanism, it is important to gain better 
control of the uncertainties of the PDFs. 
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1. Introduction 



The search for the Higgs boson(s) has been one of the central problems in high energy 
physics. It is not only important to find the Higgs boson(s), it is just as important to 
understand its properties, including couplings. These will reveal the underlying electroweak 
symmetry breaking mechanism, be it Standard Model (SM) or beyond the SM. In order 
to distinguish between different physical mechanisms, it is essential to assess the inherent 
uncertainties of theoretical predictions based on various underlying physics scenarios. 

In this regard, one of the first physical quantities that we need to gain a good theo- 
retical control of is the production cross section of the Higgs particles. A lot of work has 
been done on higher order QCD corrections to the Higgs boson production cross section 
in SM and in Supersymmetry models. These have considerably reduced the "theoretical 
uncertainties" of the calculated cross sections, usually estimated by varying the renormal- 
ization and factorization scales over some range (say, by a factor of 2). We shall refer to 
this type of uncertainties, in short, as scale uncertainties. However, the total uncertainty 
of the predicted Higgs production cross sections also includes uncertainties due to parton 
distribution functions (PDFs) — in short, PDF uncertainties. These can be significant, or 
even dominant, compared to the scale uncertainties, depending on the theoretical model 
and the model parameters. 

In this paper we perform a detailed study of the PDF uncertainty for inclusive Higgs 
boson production at the Fermilab Tevatron and CERN Large Hadron Collider (LHC) and 
compare them with existing estimates of the scale uncertainties. Particular attention is 
given to the 6-quark initiated Higgs boson production mechanism. Whereas the bb — > 
TC process is relatively insignificant compared to gg TC (via a top quark loop) in the 
SM,'^ this is not so in many models beyond the SM, such as the Minimal Supersymmetric 
Standard Model (MSSM). The presence of several vacuum expectation values in these 
models can lead to a large enhancement of the Yukawa coupling of the 6-quark to some of 
the Higgs bosons. These scenarios offer attractive opportunities for Higgs boson searches 
at hadron colliders. It is therefore important to have reliable estimates of all uncertainties 
in the calculation of the 66 — > W cross section. 



In Section 2.1, we describe general features of the models with enhanced bbTi coupling. 



as well as the specific model that is used in our numerical study — the production of the CP- 



odd Higgs boson A in MSSM. In Section |2.2| , we discuss the relevant (scheme-dependent) 



QCD subprocess for 66'H-coupling-enhanced Higgs production in these models; and the 
particular scheme used in our PDF uncertainty study. In Section |^, we present the results 
on these uncertainties for the inclusive Higgs boson production at Run II of the Tevatron 
and the LHC, and compare them to the scale uncertainties available in the literature. We 
find that the PDF uncertainties dominate over the scale uncertainties, except for low Higgs 
masses at LHC. The PDF uncertainties are calculated by the robust Lagrange multiplier 
method within the CTEQ global analysis framework, and compared to those obtained by 



^In models with more than one neutral Higgs boson, the symbol Tl is used in the generic sense to 
represent any one of the neutral Higgses, such as {h, H, A\ in MSSM, cf. Sec. 2. 
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the (more approximate) Hessian method used before. We show that the results of these 
two methods are consistent with each other. Section ^ states our conclusions. 



2. Higgs Production in the SM and in MSSM 



2.1 bbH Yukawa Coupling Enhanced Higgs Production 




H 
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Figure 1: Higgs production by: (a) bb annihilation, and (b) gg fusion through a top-quark loop. 

In the Standard Model, Higgs production by the partonic process bb ^ H is small com- 
pared to production hy gg H , cf. Fig. |^, because the bbH coupling is small compared 
to the ttH coupling, and because the 6-quark distribution is much smaller than the gluon 
distribution. However, the two production mechanisms can become comparable in beyond 
the Standard Model theories in which the bbH. Yukawa coupling is substantially enhanced 
with respect to ttTi. [1-3]. For example, in the MSSM, there are two Higgs doublet super- 
fields, with two independent vacuum expectation values (VEVs): and Vd- While the sum 
of the squares of these VEVs is fixed by the well-known Z boson mass, their ratio, denoted 
by tan f3 = Vu/vd, is a free parameter of the model. After spontaneous symmetry breaking, 
there are 5 physical particles in the Higgs sector: h (light), H (heavy), A (pseudoscalar) 
and i?^ (charged) . An important feature of this model is that with high values of tan /3, the 
Yukawa couplings of the 6-quark to the neutral Higgses are enhanced by a factor 1 / cos /3 
compared to their SM value. The t- and b- quark Yukawa couplings can be written as 
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(2.1) 
(2.2) 



and the MSSM Yukawa coupling of the Higgs bosons to t- and b- quarks, Y-^jf, relative to 
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the SM one Y^^j = gruf /2mw = mj/v, takes the form 
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(2.4) 



where a is the mixing angle of two CP-even Higgs bosons, and the weak scale v equal to 
246 GeV. 

In the MSSM, the Higgs mass is a function of niA and tan /?. The relatively high lower 
limit on the Higgs boson mass deduced from LEP data [4] favors scenarios with high tan /3. 
Also, theoretically, high tan/3 scenarios are highly motivated by SO(IO) SUSY GUTs (see 
e.g. Refs. [5-22]). Thus, it is important to explore the phenomenological consequences of 
Higgs production by enhanced bbTC Yukawa coupling mechanisms at the Tevatron and the 
LHC. 

To make our study less dependent on SUSY parameters, we focus on production of 
the CP-odd Higgs particle A. As seen from Eq. ( |2.3| & ^), the Yukawa couplings of A to 
the heavy quarks are independent of the Higgs mixing angle a; and the Abb (Att) coupling 
is enhanced (suppressed) by the factor tan /3. For simplicity, we shall not consider SUSY- 
QCD and SUSY-EW corrections to the Yukawa couplings, which could be significant in 
specific regions of SUSY parameter space. ^ 

2.2 Production Mechanisms with Enhanced bbH Coupling 

The simplest partonic processes contributing to inclusive Higgs production with enhanced 
bbH coupling are represented by the tree diagrams of Fig. (a) bb H; (b) gb Hb; 
and (c) gg — > bbTC. In QCD, these are not independent production mechanisms, since 
6-partons inside the hadron beam/target arise from QCD evolution (splitting) of gluons, 
and gluons radiate off quarks [29-31]. The three processes (a,b,c) all give rise to the same 
hadronic final states, with two B-mesons appearing in different, but overlapping, regions 
of phase space — either as beam/target remnants or as high pT particles. The distinction 
between the three processes depends very much on the factorization scheme adopted for 
the QCD calculation."^ 

For example, in the (fixed) 4-flavor scheme which is often used in 6-quark production 
calculations, there are no 6-partons by definition; hence (a) and (b) are absent, so the gluon- 
fusion process (c) is the only one contributing. By comparison, in the 5-flavor scheme, all 

''One should notice that squark contributions to gg could be important [23-26]. Even the gg 
process under study (which does not receive squark contribution at leading order if CP is conserved) can 
have sizable SUSY-QCD and SUSY-EW corrections to the Yukawa couplings, which could enhance both 
gg and bh ^ A processes (see e.g. Refs. [27, 28].) 

'^For a detailed explanation and complete references, see the appendix of [32]. 
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Figure 2: Higgs production mechanisms with enhanced bbH coupling: tree diagrams. 

three subprocesses contribute; and, in fact, all three are numerically comparable in spite 
of the differences in the nominal number of powers of ag present (0,1,2 respectively).'^ 
This is because the magnitude of the gluon distribution g{x, Q) is much larger than the 
heavy quark distribution b{x, Q) — by at least a factor of at currently available energy 
scales — hence compensating for the factors. In this scheme, when all three subprocesses 
are included, the overlap between them needs to be properly taken into account. This is 
most conveniently done by taking the m;, — > limit in the QCD calculations, while keeping 
7^ only for Yukawa couplings. 

These Hhb processes have been extensively studied in recent years for MSSM and for 
other beyond SM scenarios with similar enhanced WH couplings [33-45]. Calculations 
in the 5-fiavor scheme have been carried out to the 2-loop level [41], which considerably 
reduced the theoretical uncertainty due to the perturbative expansion, as estimated by the 
residual scale dependence. Comparison of results obtained in the 4- and 5-flavor schemes 
has also been carried out [45] . It shows consistency between the two schemes in the energy 
region of the Tevatron and the LHC. 

For the purpose of this paper — assessing the range of uncertainties associated with 
input parton distribution functions — it suffices to calculate the production cross section 




Figure 3: NLO partonic subprocesses included in the calculation of bhA coupling-enhanced Higgs 
production cross sections. 



"^This conclusion holds, except at asymptotic energies, beyond that of foreseeable accelerators, when the 
heavy quark parton distributions become comparable to that of the light quarks. 
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of the Higgs particle A, p p A + X, to the one-loop level in QCD in the 5-flavor 
scheme, which will be referred as next-to-leading order (NLO) cross section in this work. 
The Feynman diagrams representing the partonic subprocesses included in the calculation 
are shown in Fig. ^. The numerical calculation has been carried out with the program 
developed in [46], where the QCD-improved (running) Yukawa couplings have been used. 

2.3 Cross Sections for Higgs Production in SM and MSSM 

To make the discussion concrete, we now compare the order of magnitudes of Higgs boson 
production cross sections for bb annihilation and gg fusion through a top loop, gg h{A), 
in the SM and MSSM. The bb — > h{A) process is calculated as described above. The process 
gg h{A) is calculated using the HIGLU program [47], which includes the diagrams of 
Fig. |. 




Figure 4: Representative Feynman diagrams for the partonic processes gg ^ A through a top- 
quark loop that are used in the HIGLU calculation [47]. 

The cross sections for the Tevatron and LHC as a function of the Higgs mass are shown 
in Fig. ^ and Fig. ^, respectively. In these figures, the dashed lines represent the cross 



(a) Tevatron, Vs = 1.96 TeV, SM (b) Tevatron, Vs = 1.96 TeV, MSSM, fanp=rO 
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Figure 5: NLO cross sections at the Tevatron for: (a) bb ^ H and gg ^ H processes as well as 
their sum in SM; and (b) bb A and gg ^ A , Supersymmetric CP-odd Higgs boson production 
in MSSM with tan/3 = 10. 



-6- 



sections for the bb — > H{A) process; the dotted hnes the gg H{A) process; and the 
sohd hnes the combined results. Both renormaUzation and factorization scales are set to 
be Mh(a), and the PDFs used are the CTEQ6M set [48]. For the MSSM case, the results 
correspond to tan/3 = 10 and mtop = 178 GeV. 



(a) LHC,^s=14 TeV,SM (b) LHC,^s = 14 TeV, MSSM, tan^=10 
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Figure 6: NLO cross sections at the LHC for: (a) bb ^ H and gg ^ H processes as well as their 
sum in SM; and (b) bb ^ A and gg ^ A , Supersymmetric CP-odd Higgs boson production in 
MSSM with tan/3 = 10. 

One can see that in the SM case, the contribution from bb ^ H is negligible compared 
to gg ^ H . In contrast, the contribution from the Supersymmetric bb ^ A process 
becomes important even for the moderate value of tan (3 ~ 10. Except for the low Higgs 
mass region Mh < 115 GeV, the bb ^ A process is the dominant production mechanism. 
The ratio of bb ^ A to gg A processes is qualitatively similar at higher values of tan /3 
while the absolute value of the cross sections scales as tan^ (3. Relative ratios of the bb ^ A 
and gg A processes at the Tevatron and LHC are very similar, while the absolute values 
of the production rate at the LHC is about two orders of magnitude higher then those at the 
Tevatron. The cross section for both processes is enhanced with high values of tan/3 and 
can be really large. Therefore, these processes could be useful for precision measurement of 
Yukawa couplings as well as bottom-quark distributions. This underlines the importance 
of understanding the PDF uncertainties of the cross sections for these processes. 

3. Uncertainties of the Cross Sections for Higgs Production 

The PDF uncertainties for Higgs production are most reliably assessed by the Lagrange 
multiplier (LM) method [49,50]. It incorporates the calculated values of the Higgs cross 
section a in the global analysis, using the classic Lagrange multiplier technique. This 
determines the full allowed range of variation of a over the PDF parameter space. Our 
main results are obtained by this method. An alternative approach that can be applied 
without performing dedicated global analysis is the Hessian method [49,51], which has 
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been applied to the Higgs production gg ^ H process in [52] and to gg Hbb in [53]. 
The Hessian method is less robust because it relies on a linear approximation in the error 
analysis. As a part of our study, we will compare results obtained by the two methods. A 
brief summary of both methods is given in the Appendix. 

3.1 PDF Uncertainties for the 66 — > A process: the Lagrange Multiplier Method 

For a given Higgs mass, the result from a LM study of the range of variation of the Higgs 
cross section is presented in Fig. ^. The plot shows the goodness-of-fit of the global analysis 
(as measured by an overall effective value) in constrained fits, as a function of the Higgs 
cross section (for the bb ^ A process) over a certain range around the best-fit value. The 
curves represent smooth interpolations of these constrained fits (cf. the Appendix, around 
Eqs. (4.5-4^), for explanation of the method). Fig. |^(left) presents results for Tevatron for 



Ma=WO GeV, and Fig. ^(right) presents results for LHC for M4=400 GeV. These curves 
are quite close to being parabolic. This suggests that the alternative Hessian method may 
be a reasonable approximation. 
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Figure 7: Estimating the uncertainty of the Higgs cross section by the Lagrange Multiplier 
method: of the global analysis versus (j{bb — > A). The left plot presents results for Tevatron; 
the right plot for LHC. 

The uncertainty range of the Higgs cross section a is obtained by adopting a reasonable 
tolerance for the global x^, = Ax^- Various global analysis groups (CTEQ, MRST, 
ZEUS, HI) have adopted values of in the range 50-100, for ~ 2000 data points. We shall 
choose the more conservative value = 100, which we interpret as a 90% CL uncertainty 
range [48]. We obtain a+ and C7_ as the two solutions of the equation x^[o"] = T; and 
define the PDF uncertainty of the total cross section as 

Aa = {a+ - a-)/2, (3.1) 

and the relative uncertainty as 

6a = Aa/a. (3.2) 
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We evaluated this relative uncertainty of the Higgs cross section, due to the input 
PDFs, for various values of the Higgs mass. The results are presented in Fig. ^(left) 
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Figure 8: Results on the relative PDF uncertainty (in percentages) <5<t^^^ for the Tevatron (left) 
and the LHC (right) as a function of Higgs boson mass. For comparison, the corresponding QCD 
scale mrcertainty S^^^o (fj-Qj^^ j^gf [4^]^ represented by dashed lines. 

for the Tevatron, and Fig. ^(right) for the LHC, as solid lines. For comparison, QCD 
scale uncertainties available for the Higgs mass range below 300 GeV from Ref. [41] are 
represented by dashed lines. 

As one can see, there is a qualitative difference in the behavior of 6a^^^ as a function 
of Higgs mass between the Tevatron and LHC results: at the Tevatron, 6a^^^ always 
increases with increasing Higgs mass; while at the LHC, it has a minimum for Higgs mass 
around 300 GeV. To understand the reason for this behavior one can look at the uncertainty 
of the gg luminosity function, which is directly related to the 6-quark PDF uncertainty, 
since gluon splitting creates the 6-quark parton density. This is shown in Fig. ^, for the 
Tevatron on the left and LHC on the right as a function the gluon-gluon invariant mass. 
One can see that for Ma > 100 GeV, the uncertainty of gg luminosity always goes up with 
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Figure 9: Fractional uncertainty of the gluon-gluon luminosity functions at the Tevatron and 
LHC as a function of gluon-gluon invariant mass v^- 
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the increasing Higgs boson mass, which 

—V 35 I 

is related to the fact that at the Tevatron ^ 

the required x-value of PDF is already as big ^ 30 i Tevatron 

as 0.05 {x = Ma/VS > 100/1960). There- - - - lhc 

fore with the increasing Higgs boson mass, : 
X goes up and so does 6a^^^ . At the LHC, 
however, for low Ma — 100 GeV, Xmin — 
0.007 and therefore 6a^^^ is still big since 
X is fairly small. When the Higgs mass in- 
creases and reaches Ma — 300 GeV, 6a^^^ 
takes the ~ 0.02. With 

further Higgs mass increase, 6a^^^ grows ~'~ 

similarly to its behavior at the Tevatron. 

10 10 

Actually, the x-value is the principal vari- x=Mf/CME 
able that controls the PDF uncertainty. This 

can be clearly seen from Fig. which presents Figure 10: Relative PDF uncertainties ^cr^^^ 
^fjPDF ^ function of x for Tevatron and of the NLO bb A production cross sections, as 
LHC: in the x-region where Tevatron and a function of x for Tevatron and LHC. 
LHC overlap, their PDF uncertainties are 
in good agreement. 

3.2 Comparison of PDF to Scale Uncertainties 

Included in Fig. |8| is a direct comparison of PDF uncertainties 6a^^^ of the Higgs cross 
section with the scale uncertainty 6a^'~' . The latter has been obtained from the QCD scale 
dependence of the next-to-next-leading order (NNLO) calculation of the inclusive bb ^ A 
+X processes [41]. It was found that the scale uncertainty goes down from 15% to 5% at 
the LHC and from 10% to 3% at the Tevatron when the Higgs mass increases from 120 
to 300 GeV. We notice the opposite trend of those uncertainties versus Higgs boson mass 
at the Tevatron: 5a^^^ goes up from 11% to about 30% for Ma increasing from 100 to 
200 GeV, while 6a^^' decreases from 11% to 6%. Therefore, at high Ma values, 6a^^^ 
becomes almost an order of magnitude larger than ^o"^'^! 

At the LHC, both 5a^^^ and 6a^'-^ decrease with the increasing Ma in this mass 
range (100 — 200 GeV); but da^'-' is larger than 5a^^^ by a factor 1.5 — 3, depending on 
the Higgs mass. This plot suggests that at higher values of Ma, say > 300 GeV, the PDF 
uncertainties will become dominant, similar to the situation at the Tevatron. However, 
NNLO scale uncertainties were not published for this Higgs mass range. 

3.3 Comparison of the Lagrange Multiplier and Hessian Methods for Estimat- 
ing Uncertainties 

We calculated the PDF uncertainties of the cross section using the LM method, which is 
the most reliable one available. But it requires the full machinery of global analysis. The 
alternative Hessian method, utilizing a general set of eigenvector PDF sets that embody 
the PDF uncertainties [48], is more approximate, but more convenient [49,51]. Fig. pT| 
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Figure 11: Relative PDF uncertainty of tlie bb A process obtained by the Hessian method 
(dash line), compared to that by the LM method (solid line), for Tevatron(left) and LHC (right). 
Scale uncertainties from Ref. [41] are presented by the dot-dash lines. 



presents the results on PDF uncertainties of the bb ^ A cross section obtained by the 
Hessian method, compared to that obtained by the LM method for Tevatron (left) and 
LHC (right). As one can see, the two results are in good agreement. In Fig. |l^ we present 
results analogous to Fig. 11 but for the gg A process. Again, one can see that there is 
good agreement between the two methods. 
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Figure 12: Relative PDF uncertainty of the gg ^ A process obtained by the Hessian method 
(dash line), compared to that by the LM method (solid line), for Tevatron (left) and LHC (right). 
Scale uncertainties from Ref. [41] are presented by the dot-dash lines 
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3.4 Comparison of PDF Uncertainties of the bb ^ A and gg A Processes 

Notice that PDF uncertainties for the gg A process in comparison with bb ^ A are 
about a factor of two smaller. To understand this fact we refer the reader to [54] where 
Sg/g and 6b/b uncertainties and their correlation were studied in detail. For the Tevatron, 
PDF uncertainties vary from 5 to 15 — 20% for ranging between 100 and 200 GeV, and 
dominate the scale uncertainty only for heavy Higgs of mass about 150 — 160 GeV. For 
the LHC, PDF uncertainties are 5 — 6% at Ma = 100 GeV, decreasing to the minimum 
of 3 - 4% at Ma = 300 GeV and increasing again up to about 11% at Ma = 1000 GeV. 
Available scale uncertainties for Ma < 300 GeV are about a factor of two bigger than PDF 
uncertainties. Our results on gg A PDF uncertainties are in agreement with results 
presented in Refs. [52, 53]. 

We note that the PDF uncertainties for bb 
A and for gg A are strongly correlated, as 
shown in Fig. both for the Tevatron and for 
the LHC. This is hardly surprising, given the fact 
that the 6-quarks are radiatively generated from 
the gluon in the way all current parton distribu- 
tion functions are calculated. 
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Figure 13: Correlation of PDF uncer- 
tainties for bb ^ A and gg A processes 
for Tevatron and LHC. 



Finally, Fig. 14 illustrates the NLO cross sec- 
tion bands for PDF uncertainties, shown in green, 
overlaid with the NNLO scale uncertainties, shown 
in red. One can see that for low Higgs mass, 
the scale uncertainties are comparable or domi- 
nant for both colliders and both processes, while 
for heavier Higgses, PDF uncertainties are signif- 
icantly larger than the scale uncertainties. 

4. Conclusions 

The role oi bb ^ A and gg A processes may be central for the Higgs boson search. 
Therefore the correct understanding of uncertainties of their production rate is crucial. 

We found that the PDF uncertainty of 66 — > A is about a factor two larger than 
the PDF uncertainty for gg A . It was found that at the Tevatron, PDF uncertainty 
dominates the scale uncertainty for Ma >130 GeV and could be as large as 30% for 
Ma=200 GeV, which is an order of magnitude larger than the NNLO scale uncertainty. At 
the LHC the scale uncertainty is dominant and could be as big as 15% for Ma <300 GeV. 
In this region one could expect large Higgs production rates that would statistically allow 
the precision measurement of Higgs Yukawa couplings. Therefore, higher order corrections 
would be necessary in this case for better theoretical control of the cross section. For 
Mh >300 GeV, PDF uncertainty is likely to dominate at the LHC, similarly to the picture 
for the Tevatron. These results underline the importance of gaining better control of the 
PDF uncertainties, in the study of Higgs physics in the next generation of Colliders. 

We have also found that the Lagrange Multiplier and Hessian methods for assessing 
PDF uncertainties are in good agreement with each other. 
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Figure 14: NLO cross section bands for PDF uncertainties, shown in green, overlaid on NNLO 
scale uncertainties, shown in red. The upper two frames show results for bb ^ A a,t Tevatron (left) 
and LHC (right), while lower two frames show results for gg ^ A process at Tevatron (left) and 
LHC (right), respectively. 
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APPENDIX 

In this appendix we briefly review the methods of global PDF fit and estimation of PDF 
uncertainties using Hessian and Lagrange multiplier methods. 

Parton distributions which are being used for the SM and new physics predictions are 
obtained from global analysis using a "best-fit" paradigm for which the PDF is selected for 
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the minimum of the chosen ^ function. The main question is what are the uncertainties 
of those PDFs? 

In our study we are using two methods for this purpose, namely Hessian method [49, 51] 
and method of Lagrange multipUer [50] which as it was discussed in Ref . [48] , overcomes 
various problems of standard error analysis. In particular, in Ref. [48], the authors pre- 
sented a reliable way of understanding the behavior of function in the neighborhood of 
the global minimum, providing the way of correct understanding of the PDF uncertainties 
in the prediction of the cross sections. 

We summarize here briefly both methods. Both methods use a chi-square function 
is defined by 

2 2/ \ 1 2/ \ [-Ci — X^fc — Tj(a)]^ 2 i a -i\ 

X =2^Xe(a>0> where Xe(a>0 = 2^ ' 2 + 2^''fc>(4-l) 



a,- 

where e labels an experimental data set and i labels a data point in each particular data 
set. Di is the data value, aj is the uncorrelated error, and is the /cth correlated 
systematic error; these numbers are published by the experimental collaboration. Ti{a) is 
the theoretical value, a function of a set of n PDF parameters, {ai, . . . ,0^}. Also, {rfe} 
is a set of Gaussian random variables and rfc/3fcj is a (correlated) shift applied to Di to 
represent the fcth systematic error. We minimize the function ^{a^ r) with respect to both 
the PDF parameters {a} and the systematic shift variables {r^}. The result yields both 
the standard PDF model with parameters {ao}, and the optimal shifts {ffc} to bring theory 
and data into agreement. This minimum of x^ represents the best fit to the data [48]. 

Hessian method for analysis of PDF uncertainty in the neighborhood of the minima 
of involves the Hessian matrix 

_ 1 a'xo 

calculated at the minimum of Xo- The next step is to diagonalize Rij and to find its 
eigenvectors. Then for each eigenvector we have two displacements from {ao} (in the + 
and — directions along the vector) denoted {a^} and {a~} for the ith eigenvector. At 
these points, x± = Xo + ^ parametrizes the tolerance. The appropriate choice of 

tolerance T cannot be decided without a further, more detailed, analysis of the quality of 
the global fits. After studying a number of examples [50, 51], we concluded that a rather 
large tolerance, T ~ 10, represents a realistic estimate of the PDF uncertainty. 

One can show that in a linear approximation, the uncertainty 5X for any quantity X, 
which depends on PDF, can be expressed as 

,,Xf ,4.3) 
or, in terms of the eigenvector basis sets. 



1 " 

{bXf = -Y^{X{at)-X{a-)\\ (4.4) 



4 
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representing the master equation defined in [50]. One should point out again, that equation 
( |4.4| ) is based on a hnear approximation: x^(^) is assumed to be a quadratic function of 
the parameters {a}, and X{a) is assumed to be linear. This approximation is not strictly 
valid in general. 

The essence of the Lagrange multiplier method is the introduction of the Lagrange 
multiplier variable A and minimizing the function 

xl{X,a)=x\a) + XX{a) (4.5) 

with respect to the original n parameters {a} for fixed values of A. In Eq j4.5| X(a) is some 
observable as in the example for Hessian method. Minimization of Xxi^ju) for various 
values of A allows to find the parametric relationship between x^(o^) a-^d X{a), .i.e. 

x!(A, ao) = xHao) + XX{ao) ^X = X{x\ao, A)), (4.6) 

where oq is the set of parameter values {a} for each particular value of A. Eq. ( [4.6[ ) is the 
key point of LM method since for given value of tolerance 

Ax^ = x\ao,X±)-x\ao,0) (4.7) 

which would correspond to some two values A^, one can find the respective variation of 
the observable X: 

6X+ = X{x\ao,X^)) - X{x\ao,0)), 5X^ = X{x\ao, X^)) - X{x\ao,0)) . (4.8) 

The LM method for calculating 5X± is more robust in general since it does not approxi- 
mate X{a) and x^(fl) by linear and quadratic dependence on {a}, respectively, around the 
minimum. 
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